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Abstract 

The physical shape of a giant planet can reveal important information about its centrifugal potential, 
and therefore, its rotation. In this paper I investigate the response of Jupiter's shape to differential 
rotation on cylinders of various cylindrical radii using a simple equipotential theory. I find that both 
solid-body rotation (with System III rotation rate) and differential rotation on cylinders up to a latitude 
between 20 and 30 degrees are consistent with Jupiter's measured shape. Occultation measurements of 
Jupiter's shape could provide an independent method to constrain the depth of its zonal winds. 

1 Introduction 

Despite decades of investigation, Jupiter's atmospheric circulation is still not well understood. Major ques- 
tions revolve around Jupiter's dynamical structure, in particular, how its zonal winds are forced (from below 
by internal heat or from above by solar heating) and how deeply they extend into the interior. Many au- 
thors have attempted to better constrain and understand Jupiter's dynamics based on various approaches 
(typically deep convection models vs. shallow atmospheric models) and different numerical methods (e.g., 
Williams, 1978; Cho and Polvani, 1996; Nozawa and Yoden, 1997; Huang and Robinson, 1998; Williams, 
2003; Aurnou and Heimpel, 2004; Heimpel et al., 2005; Vasavada and Showman, 2005; Showman et al., 2006; 
Heimpcl and Aurnou, 2007; Lian and Showman, 2008). 

Deep convection models of the Jovian zonal winds (e.g., Busse, 1976; Heimpel et al., 2005; Heimpel and 
Aurnou, 2007) suggest that the bottom boundary of the deep convection is at about 0.9 Rj (Rj is the radius 
of Jupiter) . Liu et al. (2008) discuss the importance of Lorentz forces in constraining the depth of the Jovian 
circulation. They argue that the observed zonal winds on Jupiter cannot extend more deeply than 0.96 Rj. 
Shallow forcing models typically argue that the zonal winds are driven by moist convection and are a surface 
phenomenon, i.e., a 'weather layer' (e.g., Ingersoll et al., 1969; Ingersoll et al., 2004). 

The penetration depth of Jupiter's zonal winds also has implications for the planet's internal structure. If 
the zonal winds are 'deep' and involve a non-negligible amount of mass, interior models that typically assume 
solid-body (SB) rotation must be modified, and corrections to the gravitational coefficients due to dynamics 
must be included (Hubbard, 1982; 1999). A recent interior model of Jupiter by Militzer et al. (2008) suggests 
that Jupiter's measured J4 value can only be fit when differential rotation is considered because no interior 
models that fit Jupiter's J4 could be found under the assumption of solid-body rotation. The interior model 
presented by Militzer et al. (2008), however, consists of 2-layers (a heavy-element core surrounded by a 
hydrogen- helium envelope) , while standard 3-layer Jupiter models with solid-sbody rotation can typically fit 
Jupiter's gravitational coefficient (J2,J4,Ja) (e.g., Saumon & Guillot, 2004; Nettelmann et al., 2008). As 
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the question of whether Jupiter's interior structure is consistent with solid-body rotation remains open, it is 
clear that a better determination of the depth of Jupiter's zonal winds is desirable for both dynamical and 
interior models. 

Jupiter's internal dynamics, or more precisely, its rotation profile, can be constrained by measurements 
of its high-order gravitational coefficients. Hubbard and collaborators (e.g., Hubbard, 1982; 1999; Kaspi, 
et al., 2010) showed that a departure from solid-body rotation is detectable in gravitational coefficients 
larger than about degree 10 (Kaspi et al., 2010, figure 7b). The Juno mission to Jupiter (Bolton, 2005) is 
designed to provide accurate measurements of Jupiter's gravitational field to high order and can therefore 
constrain its dynamical structure. Gravitational anomalies can also be used directly to infer Jupiter's density 
anomalies and therefore its rotation profile (see Kaspi et al., 2010 and references therein). In this paper I 
use an equipotential theory to derive Jupiter's shape when differential rotation on cylinders is included 
and suggest that shape data, i.e., occultation radii, can provide an independent method of constraining its 
internal rotation profile. 

2 Jupiter's Geoid with Differential Rotation on Cylinders 

The effective potential U of a giant planet in hydrostatic equilibrium is given by a combination of its 
gravitational potential and its centrifugal potential (Kaula, 1968), 

U = — ~ g (r) J2 " P2 " (M) J + 2 q {— ) (a) (1 - (1) 

where GM is the gravitational constant times the total planetary mass, a is the equatorial radius, Ji n are the 
gravitational coefficients, and P^n are the Legendre polynomials. The point where the potential is evaluated 
is given by the planetocentric radius r and latitude 4> {p- = smcf>). The second term on the right side is the 
centrifugal potential, q = uj 2 a 3 /GM is the smallness parameter, where to is the angular velocity of rotation. 
The planet can rotate uniformly as a solid-body, or on cylinders to a certain depth, including differential 
rotation on cylinders all the way to the planet's center. The harmonic coefficients J 2n are obtained from the 
measured values for a reference equatorial radius of 71,492 km (Jacobson, 2003). The values of Jm from 
Jacobson (2003) in units of 10~ 6 are J 2 = 14696.43 ±0.21, J 4 = -587.14 ± 1.68, J 6 = 34.25 ± 5.22. 

Helled et al. (2009) discuss how the planetary shape derived from the planet's effective potential, can 
be used to gain information on its internal rotation. Helled et al. (2009) searched for a solid-body rotation 
period that best fits Jupiter's measured physical shape, i.e., a solid-body rotation period that minimizes 
the planet's dynamical heights (for details, see Helled et al. 2009 and Lindal, 1992). The physical shape of 
Jupiter was derived by adding the dynamical heights (obtained using wind data) to the calculated geoid. The 
shape of Jupiter's geoid in Helled et al. (2009) was derived assuming solid-body rotation, including a range of 
rotation rates. Since Jupiter may not rotate solely as a solid-body, it is desirable to perform an investigation 
of Jupiter's geoid shape in which the possibility of differential rotation is included. In the analysis below 
I project Jupiter's zonal winds along cylinders to describe its rotation profile, which is explicitly used in 
the reference effective potential equation (Eq. 1). The derived shape using different cylindrical radii is then 
compared with Jupiter's measured shape. Figure 1 shows Jupiter's physical shape as calculated by Helled et 
al. (2009) with six measured occultation radii from the Pioneer and Voyager missions (see Fig. 7 in Lindal 
et al., (1981). 



2 



When differential rotation is included, the planetary rotation rate is no longer given by a constant (u>o); 
instead, the rotation profile changes with cylindrical radius. I assume that Jupiter rotates on cylinders about 
the polar axis, so the angular velocity profile ui depends on the distance of a cylindrical shell from the spin 
axis (Helled et al., 2010). I define 7 as the normalized cylindrical polar radius, and derive the angular velocity 
w(7) using Jupiter's wind data. I use Jupiter's zonal wind velocities obtained from Hubble Space Telescope 
images from 1995 to 2000 (Garcfa-Mclendo & Sanchez-Lavega, 2001) at three different wavelengths (410, 
892, 953 nm) representing Jupiter's cloud top levels (~ 1 bar). The value of 7 ranges between zero and one; 
7 = 1 corresponds to solid-body rotation, and 7 = refers to differential rotation on cylinders all the way to 
Jupiter's center. Figure 2 shows a sketch of the planet's rotation profile for a given normalized cylindrical 
polar radius value (70). The planet rotates as a solid-body (with a rotation rate uj n ) up to 70, from that 
point on the planet rotates differentially on cylinders with a rotation rate profile 0^(7). 4>^ diSS is the latitude 
corresponding to 70 that defines the latitude in which transition between solid-body and differential rotation 
occurs. 

Jupiter's zonal wind velocities are available for both the northern and southern hemispheres. Although 
the wind profiles are similar in the two hemispheres, they are not identical. Figure 3 shows uj(-f) using the 
northern (dashed black curve) and southern (dashed gray curve) wind velocities separately. The horizontal 
thin black line shows the value of ujq, Jupiter's angular velocity for a solid-body rotation period of 9h 55m 
29.56s (Higgins et al., 1996). Note that 7 = corresponds to the planet's pole (cj) UJdiff =90 o , i.e., differential 
rotation on cylinders all the way to the center) while the wind velocity data do not reach the poles. In 
the northern hemisphere, wind velocities are measured to higher latitudes and therefore can reach smaller 7 
values. The thick black curve shows the angular velocity averaged over the northern (dashed black curve) and 
southern (dashed gray curve) hemispheres. I can now compute Jupiter's centrifugal potential for different 
values of 7, i.e., differential rotation on cylinders for different cylindrical polar radii. 

Reference geoids representing Jupiter's shape (Helled et al., 2009), are defined in the usual manner as the 
level surfaces of equal effective potential, and are derived for both the 100 mbar and 1 bar pressure-levels. 
The advantage in using the 100 mbar level is that radio occultation data refer directly to this pressure- 
level. The zonal wind velocities, however, correspond to ~ 1 bar. As shown below, I get similar results 
for both pressure-levels. Jupiter's polar radius is held fixed to its measured value and is used to define the 
cquipotential surface. The polar radii for the 100 mbar and 1 bar pressure levels are 66,896±4 km and 
66,854±4 km, respectively (Lindal et al., 1981; Lindal, 1992). The geoid's radius is calculated to fifth order; 
the harmonics J2, J4, and Jq are extrapolated to even harmonics Jg and J 10 to ensure that the precision 
of the geoid calculation is better than ± 1.0 km. There is an implicit assumption in the procedure that 
the mass in the atmosphere above the calculated pressure level is negligible. The computation therefore 
assumes that non-hydrostatic effects have a negligible contribution to the planetary shape. More details on 
the method can be found in Helled et al. (2009, 2010) and references therein. 

3 Results 

The six radio occultation radii of Jupiter's 100 mbar isosurface at different planetocentric latitudes, four of 
them at the northern hemisphere, are shown in Fig. 1. The formal error on Jupiter's measured radii is ^4 km 
(Lindal et al, 1981; Lindal, 1992), and can therefore be taken as the standard error bar on Jupiter's shape. 
Helled et al. (2009) presented the extrapolated occultation radii at the 1 bar pressure level by subtracting 
47.82 km along the vertical to the geoid and projecting it on the geocentric radius. The estimate of the 47.82 
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km height of the 100 mbar level above the one-bar level in the atmosphere is obtained by interpolating in 
Table 8.2 in Lodders & Fegley (1998). 

Below I vary the differential rotation polar cylindrical radius and investigate how Jupiter's shape is 
changing accordingly. I compute Jupiter's shape assuming that differential rotation on cylinders occurs 
up to each available occultation measurement at a given planetocentric latitude cj>. I then compare the 
calculated shape with the measurements. The planetary shape is adjusted so the radii at the transition 
between differential and solid-body rotation converge (Helled et al., 2010). I.e., for planetocentric latitudes 
greater than <fiu> diff the planetary shape is similar to that of Jupiter with a solid-body rotation. Although 
this assumption ensures that the planetary shape is an equipotential surface, it neglects the zonal winds 
at latitudes larger than <f>w iiff . As a result, configurations in which differential rotation penetrates only 
partially at high latitudes (e.g., Hiempel et al., 2005) are not considered. The geostrophic balance on jets 
at high latitudes would force the 100 mbar (or 1 bar) pressure-level to exhibit topography relative to the 
equipotential, and therefore solid-body rotation would lead to smaller radii at these latitudes. However, the 
amplitude of this effect is ~ 1-2 km, smaller than the formal error bar on Jupiter's radii. 

3.1 Averaging North and South 

Table 1 shows the calculated radii for the occultation latitudes at both the 100 mbar (top) and 1 bar (bottom) 
pressure-levels when the zonal winds are averaged between north and south. r occ is the measured radius at a 
given pressure- level. Ar of SB rotation gives the radii residuals in km for solid-body rotation. Ar oiuidiff at 
a given planetocentric latitude is the difference (in km) between the measured and computed radii assuming 
differential rotation up to a given <j> Udiff (see figure 2). For latitudes <j>w di!J the residuals are similar to those 
computed for SB rotation. 

Since the wind velocities are averaged between the northern and southern hemispheres, results are pre- 
sented for all occultation latitudes, and for both 100 mbar and 1 bar pressure levels. Generally, better 
agreement with the measurements, regardless of the rotation profile, is found for the 100 mbar pressure- 
level. This is due to the lack of extrapolation in radius as done for the 1 bar pressure- level (Helled et al., 
2009). The residuals found here are typically larger than the ones presented by Helled et al. (2009) because 
I do not include the offset in north-south symmetry, which is found to be <~ 5.5 km for Jupiter (for details, 
see Helled et al., 2009 and references therein). For both pressure-levels, solid-body rotation leads to shapes 
that agree well with occultation data, as found by Helled et al., (2009). Differential rotation on cylinders, 
depending on its cylindrical polar radius ('depth'), could lead to similar agreement as long as differential 
rotation does not penetrate all the way to Jupiter's center. Given that the error on Jupiter's shape is ^4 
km, residuals in radius up to about 5 km can be considered as a good fit. In this case the results suggest that 
differential rotation on cylinders above $u) diSS ~ 20° can be ruled out for Jupiter. It should be noted that 
when differential rotation up to higher latitudes is considered, the discrepancy with observations occurs at 
low-latitude regions. Differential rotation at high latitudes can lead to agreement with high-latitude shape 
but not with the equatorial shape. 

In order to better understand this behavior, one could consider the following exercise. As seen in Fig. 3, 
Jupiter's jets are mostly eastward, and thus a model with jets that penetrate through the planet essentially 
has a slightly faster rotation rate than a planet in solid-body rotation at the System III rotation rate. A 
typical speed of the mid- and high-latitude jets is ~ 20 m s _1 , about 0.2% of the speed of the planet's 
rotation in inertial space. To zero-order, increasing the rotation rate by 0.2% should increase the rotational 
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flattening by the same factor. The difference in equatorial and polar radii for Jupiter is about 4600 km. 
Thus, increasing the flattening by 0.2%, would imply an increase in the equatorial radius (relative to the 
polar radius) of the order 4600 x 0.002 = 9 km. In comparison, Table 1 quotes equatorial radii that arc 
about 11 km larger when differential rotation up to high latitudes is included. 

3.2 Separating North and South 

I repeat the calculation, this time separating the wind data for the northern and southern hemispheres. The 
results are summarized in Table 2. For the northern hemisphere, the differences between the measured radii 
and the calculated ones are slightly worse than in the previous case; the results, however, remain the same. 
Again, I find that solid-body rotation provides a good fit to wind data as well as differential rotation up to 
a latitude <pu di}1 ~ 20°. Differential rotation toward Jupiter's center results in larger discrepancy with the 
shape data. 

For the southern hemisphere, only two measurements are available, one of them is at latitude of 71.8° 
S. Since Jupiter's zonal wind velocities in the southern hemisphere are not available at this latitude, only 
the occultation measurement at 10.1° S is presented. The results suggest that both solid-body rotation and 
differential rotation on cylinders up to (j> Udiff = 10° are consistent with Jupiter's measured shape, a result 
which is consistent with the previous cases. 

4 Summary 

The centrifugal potential of a giant planet influences its shape. Knowledge of Jupiter's shape can be used 
to better constrain the depth of its zonal winds. Better determination of the zonal winds' penetration depth 
is important for the understanding of Jupiter's global circulation, as well as its magnetic field generation 
(e.g., Ingersoll et al., 2004). I present a simple computation of Jupiter's geoid when differential rotation 
on cylinders is included and compare its shape to available shape data. To keep the planet in hydrostatic 
equilibrium and treat its shape as an equipotential surface, I define differential rotation on cylinders about 
the polar axis. For this rotational configuration, I find that the maximum latitude to which differential 
rotation on cylinders can occur lies between 20 and 30 degrees. 

While the work presented here assumes differential rotation on cylinders, which may be a good represen- 
tation of Jupiter's internal dynamics, other rotational profiles may be possible, such as differential rotation 
on spheres, as seems to be the case for the Sun (Elliott et al., 2000; Glatzmaier, 2008). Although further 
investigation of the topic is desirable; in particular, shape models that account for non-hydrostatic contri- 
butions, and 2-D interior models, the analysis presented here suggests that occultation shape data could be 
used to improve our understanding of Jupiter's internal dynamics. 
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Jupiter's Physical Shape at 100 mbar 
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Figure 1: Jupiter's physical shape at 100 mbar (Helled et al., 2009). The curve represents the shape that 
best fits the radio data. The discrete points correspond to the six occultation radii measured by Pioneer 10 
and 11, and Voyager 1 and 2. 
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Figure 2: A sketch of Jupiter's rotation profile, loq is the solid-body rotation rate, 0^(7) is the rotation 
rate assuming differential rotation on cylinders. 70 is the normalized cylindrical radius that defines the 
boundary between solid-body and differential rotation. 4>u, Mft is the latitude of the transition, and essentially 
determines the cylindrical radius of differential rotation as described in the text. 
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Angular Velocity a) from HST Wind Data 
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Figure 3: Angular velocity uj as a function of normalized cylindrical polar radius. The black and gray dashed 
curves correspond to the northern and southern hemispheres, respectively. The thick black curve shows the 
average northern and southern values. For 7 < 0.4 the thick black curve is identical to the angular velocity 
of the northern hemisphere since there are no wind data available for high latitudes for Jupiter's southern 
hemisphere. 
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Table 2: Jupiter's shape as calculated from the model for different cylindrical polar radii with Jupiter's zonal 
winds for the northern (top) and southern (bottom) hemispheres. The first two columns are Jupiter's plan- 
ctocentric latitude <j> and radius r obtained by occultation measurements at Jupiter's northern hemisphere. 
The other columns give the difference in radius in km between the measured radius and the calculated radius 
assuming different rotation profiles. Column three presents the radius residuals assuming solely solid-body 
(SB) rotation. The rest columns lists the residuals in radius for differential rotation on cylinders up to a 
given latitude <j>u di/f ■ 



Jupiter's Radii (100 mbar)- Northern Hemisphere 



<t> (deg) 


r occ (km) 


Ar (km) 


Ar (km) 


Ar (km) 


Ar (km) 


Ar (km) 






SB rotation 


<t>u«ff= l ° 




4^=30° 




60.3 


67886.18 


-4.45 


SB 


SB 


SB 


-4.32 


28.0 


70367.34 


-6.51 


SB 


SB 


-6.02 


-2.08 


19.8 


70896.17 


-2.11 


SB 


-2.04 


0.96 


4.97 


0.07 


71490.79 


0.58 


0.59 


5.43 


8.50 


12.59 




Jupiter's 


Radii (1 bar)- 


Northern Hemisphere 






4> (deg) 


r occ (km) 


Ar (km) 


Ar (km) 


Ar (km) 


Ar (km) 


Ar (km) 






SB rotation 




</W w =20° 




<f>w di ff= Q1 ° 


60.3 


67886.18 


1.30 


SB 


SB 


SB 


1.42 


28.0 


70367.34 


-4.46 


SB 


SB 


-3.99 


0.04 


19.8 


70896.17 


-0.83 


SB 


-0.76 


2.18 


6.29 


0.07 


71490.79 


0.81 


0.82 


5.77 


8.77 


12.96 



Jupiter's Radii (100 mbar)- Southern Hemisphere 



-10.1 71330.92 -1.56 -1.52 
Jupiter's Radii (1 bar) - Southern Hemisphere 

-10.1 71330.92 -1.02 -0.98 
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